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Unpolarized cross sections and double-helicity asymmetries of single-inclusive positive and neg- 
ative charged hadrons at midrapidity from p+p collisions at ^/s = 62.4 GeV are presented. The 
PHENIX measurement of the cross sections for 1.0 < pr < 4.5 GeV/c are consistent with pertur- 
bative QCD calculations at next-to- leading order in the strong-coupling constant, Us. Resummed 
pQCD calculations including terms with next-to-leading-log accuracy, yielding reduced theoretical 
uncertainties, also agree with the data. The double-helicity asymmetry, sensitive at leading order 
to the gluon polarization in a momentum-fraction range of 0.05 < Xgi^on ^ 0.2, is consistent with 
recent global parameterizations disfavoring large gluon polarization. 

PACS numbers: 13.85.-t,13.88.-|-e,14.20.Dli 



I. INTRODUCTION 



The comparison of cross-section predictions with data 
on single-inclusive hadron production in hadronic col- 
lisions, p+p — > h + X, is important for understand- 
ing perturbative quantum chromodynamics (pQCD). 
For hadrons produced with transverse momenta pt ^ 
Age D , the cross section factorizes into a convolution in- 
volving long-distance and short-distance components 
[^. Long-distance components include universal parton- 
distribution functions (PDFs) describing the partonic 
structure of the initial hadrons and fragmentation func- 
tions (FFs) for the final-state hadron. The short-distance 
part describes the hard scattering of partons. The long- 
distance components, PDFs and FFs, can be extracted 
from other processes, such as deep-inelastic scattering 
and hadron production in e^e~ colliders. This allows for 
a test of the short-distance part of the convolution, which 
can be estimated using pQCD. In particular, differences 
between data and predictions can indicate the impor- 
tance of neglected higher-order terms in the expansion 
or power-suppressed contributions [sj. 

Next-to-leading-order (NLO) pQCD and coUinear fac- 
torization successfully describe RHIC cross-section mea- 
surements at a center-of-mass energy {\fs)oi 200 GeV 
for midrapidity neutral pions 4, 5|, jets ^-i^d di- 

rect photo ns p , as well as forward rapidity pions and 
kaons (lol . [lll |T However, at lower ^/s, in particu- 
lar in fixed-target experiments with 20 % y/s ^ 40 
GeV, NLO pQCD calculations significantly underpredict 
hadron production, by factors of three or more [3| . The 
consistency between NLO estimations and data at low ^/s 
was improved [3, Il2l . [isj by including the resummation of 
large logarithmic corrections to the partonic cross section 
to all orders in the strong coupling as. The corrections 
are of the form Xt?^ (l — for the fcth-order term in 

the perturbative expansion. Here xt = '^PtI^T^^ where 



Pt = Pt I ~^ is the transverse momentum of the parton 
fragmenting into the observed hadron with a fraction z 
of the parton transverse momentum, and -s/l = \fx\X2s 
is the partonic center-of-mass energy where x\,xi are 
momentum fractions carried by two interacting partons. 
The corrections are especially relevant in the threshold 
regime 5;^ — ^ 1 in which the initial partons have just 
enough energy to produce a high-transverse-momentum 
parton fragmenting into the observed hadron. In this 
regime gluon Bremsstrahlung is suppressed, and these 
corrections are large tl^]- However, the addition of the 
resummed next-to-leading-log (NLL) terms to an NLO 
calculation may not provide the best means of describ- 
ing data in a given kinematic region, for example, when 
the (unknown) higher-order terms that are omitted from 
the calculation have comparable magnitude and opposite 
sign to the NLL terms. It is therefore important to test 
pQCD calculations against data in a region of intermedi- 
ate \fs, to better define the kinematic ranges over which 
pQCD calculations can be applied with confidence. 

The PHENIX-dectector data presented here for non- 
identified charged-hadron production make use of ap- 
proximately five times the acceptance and different de- 
tection techniques compared to our previous identificd- 
charged-hadron analysis [T!] . The new results allow tests 
of NLO and NLL predictions based on an independent 
measurement. In addition, the theoretical calculations 
make use of different fragmentation functions than for 
identified particle calculations. Alternatively, assuming 
the reliability of the short-distance aspects of the theory, 
the data may be used to refine our knowledge of frag- 
mentation functions. Present measurements also cover a 
greater pT range than the identified charged-hadron cross 
sections, where the measured momentum ranges for pi- 
ons, kaons, and protons are 0.3-3 GeV/c, 0.4-2 GeV/c 
and 0.5-4.5 GeV/c respectively These cross-section 
measurements of nonidentified charged hadrons (combi- 
nations of TT^, K'^ , p^) are also important as baselines for 
extracting nuclear modification factors in high-pr hadron 
production in heavy ion collisions at RHIC [IE 



16|. 
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The charged hadrons in these measurements were pro- 
duced from collisions of transversely- and longitudinally- 
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polarized proton beams, a unique capability of 
RHIC [13 • While the cross-section measurements dis- 
cussed above require averaging over the beam polariza- 
tions, sorting the hadron yields by colliding proton helic- 
ities (for longitudinal beam polarizations) provides sensi- 
tivity to the helicity PDFs [l8| . The ability to probe he- 
licity PDFs is essential for understanding the spin struc- 
ture of the proton |19j|. 

The spin of the proton originates from the spin and or- 
bital angular momenta of its quark, antiquark, and gluon 
constituents. The contribution carried by quark and an- 
tiquark spin, determined from polarized deep-inelastic 
scattering (pDIS) experiments [20l - [3g | using polarized 
leptons and polarized nucleons, is ~ 25% [13, I3ll - l33l |. 
This is surprisingly small [loj and implies that the ma- 
jority of the spin of the proton must originate from gluon 
spin and/or orbital angular momentum. 

Colliding longitudinally-polarized proton beams pro- 
vides sensitivity to the gluon-helicity distribution func- 
tion at leading order. The helicity-dependent difference 
in hadron production is defined as: 



At = 62.4 GeV, the production of final-state 
hadrons at midrapidity in a transverse momentum range 
1.5 < < 4.5 GeV/c is dominated by quark-gluon scat- 
tering 34]. This makes the asymmetry reported here, 
All{p+P + X), sensitive to the polarized gluon 

PDF AG{x) at leading order, and more sensitive to its 
sign than processes dominated by gluon-giuon scattering 
or which produce isospin-symmetric particles. For exam- 
ple, preferential fragmentation of up quarks into positive 
pions and down quarks into negative pions, combined 
with the fact that the up quark helicity PDF is positive 
and the down quark helicity PDF is negative, would lead 
to an ordering of the asymmetries of pions (charged and 
neutral) directly sensitive to the sign of the gluon helicity 
PDF. Positive AG{x) would lead to Al^^ > Af^ > All, 
whereas a negative AG{x) would imply an opposite or- 
dering. These results can be combined with data from 
polarized collider and fixed target experiments in a global 
analysis to reduce uncertainties on the gluon helicity dis- 
tribution [ssl. [36|. 



dAc 
dpT 



da^ 



dpi 



dpi 



where the superscripts -t--|- and H — refer to same and op- 
posite helicity combinations of the colliding protons [18| . 
Factorization allows this to be written as a convolution of 
the long- and short-distance terms summed over all pos- 
sible flavors for the partonic interaction a + b ^ c + X' , 
where c fragments into the detected hadron h: 



dAa 
dpT 



J dXadXbdZc Afa{Xa,tJ-f) Afb{xb, fJ-f) 



abc 



X ■ 



dpT 



{XaPa.XbPb, P^'/Zc, i^f^tJ-'f, Mr) 

(1) 



where Af{x,^f) are the polarized PDFs of the colliding 
partons carrying light-cone momentum fraction x evalu- 
ated at factorization scale fif. The fragmentation func- 
tion of scattered parton c into hadron h with fraction 
Zc of the scattered parton momentum is D'^{zc, iJ-'f) at 
fragmentation scale /i^. The helicity-dependent differ- 
ence in the cross section of the hard partonic scattering 
a+b —>■ c+X' is denoted by dAa and is calculable in per- 
turbative QCD. Cross section calculations to finite order 
in as have a dependence on factorization and renormal- 
ization scales fif and fir- 

Instead of directly measuring the helicity-dependent 
cross-section difference dAa/dpr, we extract the double- 
longitudinal spin asymmetry defined as the ratio of the 
polarized to unpolarized cross sections A^^ = dAa /da. 
Here, da is the helicity-averaged (unpolarized) cross sec- 
tion da = [dfT(+_|_) + da(^^ The ratio dAa /da has 

smaller systematic uncertainties since some of the uncer- 
tainties cancel. 



II. EXPERIMENTAL SETUP 

This analysis uses the PHENIX central arm spectrom- 
eters. Each arm has an acceptance covering a pseudora- 
pidity range I77I < 0.35 and A(f> — ^ in azimuth [stI. [38|. 
The PHENIX central magnet creates an axial magnetic 
field with / Bdl = 0.78 T-m at f in this region. 

Midrapidity charged hadrons are tracked in the drift 
chambers (DC), which are located outside the magnetic 
field with an inner radius of 2.0 m and outer radius of 
2.4 m. The tracks are reconstructed as nearly straight 
lines and yield the deflection in the axial magnetic field 
to determine the transverse momentum with a resolution 
^ = 0.007 0.009 PT (GeV/c) The first term 

in the resolution is dominated by multiple scattering in 
material before the DC, while the second arises from the 
finite angular resolution of the DC. The momentum scale 
is set by requiring the proton mass reconstructed from 
the measured momentum and time of flight to match the 
known proton mass. 

Track reconstruction also utilizes two layers of pad 
chambers (PC), which are multiwire proportional cham- 
bers with pad readout [sl]. The first layer, PCI, is lo- 
cated after the DC, with an average radius of 2.49 ni. 
PCI information is used in conjunction with the DC hits 
and vertex information to determine the polar angle for 
each charged track. The outermost layer PC3 is at an 
average radius of 4.98 m and is used for charged track 
selection by matching PC3 hit positions with tracks pro- 
jected using information from the DC and PCI and the 
measured event vertex. Matching projected tracks with 
hit positions in PCS also helps in rejecting decay back- 
grounds from primary hadrons. 

Vertex and timing information is provided by two 
beam-beam counters (BBCs) [s^ placed around the 
beam pipe. The BBCs are located 1.44 m forward and 
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backward of the nominal interaction point. Each BBC 
comprises an array of 64 phototubes fitted with 3 cm long 
quartz radiators. The phototubes detect Cerenkov radi- 
ation from charged particles traversing the quartz. The 
detectors cover a pseudorapidity range of 3.0 < |?7| < 3.9, 
and a full — 2t: in azimuth. A coincidence of hits from 
the two BBCs forms the minimum bias trigger, with tim- 
ing information providing the location of the event vertex 
along the beam line with a few cm precision. 

To eliminate the background due to photon conver- 
sion in material before the DC (primarily the beam pipe 
and DC entrance window), the analysis uses information 
from a ring imaging Cerenkov detector (RICH) [38| lo- 
cated after the DC. The RICH uses CO2 at atmospheric 
pressure as a radiator, with a momentum threshold of 
17 MeV/c for and 4.7 GeV/c for charged pious. A 
RICH veto (ensuring no RICH hits) is used to reject e"^ 
and results in an upper transverse-momentum limit of 
Pt < 4.5 GeV/c for charged hadrons in the analysis. 

Zero degree calorimeters (ZDC), which detect neutral 
particles near the beam pipe {9 < 2.5 mrad) are used 
in conjunction with the BBC to estimate the systematic 
uncertainty on the relative luminosity for the asymmetry 
measurements The BBC and ZDC are also used 

to determine the integrated luminosity measurement [4l], 

The stable spin direction of polarized protons in RHIC 
is vertical. The spin direction can be rotated into the 
longitudinal direction at the PHENIX interaction region 
using pairs of spin rotators. The polarizations of the 
beams at RHIC are measured every few hours by carbon 
target polarimeters (isj . which are normalized to an ab- 
solute measurement with a hydrogen jet target polarime- 
ter [43. 

III. CROSS SECTION 

The results presented here are the first measurements 
of the cross section of inclusive charged- hadron produc- 
tion at midrapidity in the transverse momentum range 
0.5 < Pt < 4.5 GeV/c from p+p collisions at y/s = 62.4 
GeV. The analysis techniques are similar to methods de- 
scribed in Ref. [45| and are briefiy explained in Sect. lIII Al 
The cross section results are presented and discussed in 
Sect.UHBl 



A. Cross Section Measurement 

Approximately 2.14 x 10^ BBC-triggered events corre- 
sponding to an integrated luminosity of 15.6 nb~^ from 
polarization-averaged data taken in 2006 have been 
analyzed. We calculate the midrapidity charged- hadron 
production cross section using the following formula: 

E Sa _ gBBC (PN^pt) „ 

c #3 - N^^^ dcl^pTdpTdy'^~'^'''^E^^^ ' 



where ctbbc is the y + y cross section seen by the BBC as 
measured in Ref. [40], A^bbc is the total number of BBC- 
triggered events analyzed, i?smcar is the correction factor 
for the smearing of track px owing to the momentum 
resolution of the detectors as well as multiple scattering 
of the hadron tracks, Ctrig is the correction factor for 
BBC trigger bias and is the combined correction 

factor for geometrical acceptance of the detectors and 
reconstruction efficiency. 

The reconstructed charged tracks in the transverse mo- 
mentum range 0.5 < pt ^ 4.5 GeV/c from events with 
vertices within ± 30 cm of the nominal interaction point 
are matched to projected hit positions in PC3 in az- 
imuthal (</)) and beam direction (z). Distributions of 
the matching variables (difference between the projected 
track position and actual hits on PC3, termed PC3d(/) and 
PC3dz) are fit with the combination of a signal Gaussian 
and a second Gaussian for background. Figure [1] demon- 
strates the method for a sample pt bin. The width of 
the fit to the signal distribution is used to impose a si- 
multaneous selection window of 2a for both matching 
variables. 




-0.015 -0.01 -0.005 0.005 0.01 0.015 

PC3d(|) (rad) 



FIG. 1: (color online) Difference between the track- 
extrapolated and actual hit positions on PCS in azimuth for 
a sample pr bin. The histogram represents data; the dashed 
line represents a two-Gaussian fit. 

The background to these measurements comes from 
several sources. One source is soft electrons from the 
magnet pole faces, and another is decays in flight of tt^ 
and . Not all of these electrons are rejected by the 
RICH cut. Off- vertex electrons and daughter particles 
with a perpendicular momentum kick from the decay are 
reconstructed as apparent high-p^ tracks, but with wide 
Gaussian track matching distributions. The background 
fraction in each of the px bins is determined by using 
the distributions of the matching variables. Background 
fractions, which are 2-5% for px < 2.75 GeV/c and ^ 
30% in the highest pT bin, are subsequently subtracted 
from hadron yields. 

An additional source of background is the feed-down 
background produced by weak decays of mostly A par- 
ticles close to the event vertex with apparent momenta 
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FIG. 2: (color online) Combined acceptance + detection effi- 
ciencies for separate hadron species. 



close to their true momenta and matching distributions 
peaked under the signal. Feed-down contributions to the 
detected protons and antiprotons from weak decays of 
A's and heavier hyperons are estimated using input A 
and A spectra from p + p measurements at ^/s — 63 GeV 
at the ISR \M, M and at = 62.4 GeV at PHENIX 
with a GEANT3 -based simulation of the PHENIX de- 
tector. The fractional contributions of the feed-down 
protons and antiprotons are independent of pt above 
Pt = 2 GeV/c and are close to 7 and 15% respectively. 
Below Pt = 2 GeV/c the fractions increase with decreas- 
ing Pt and are roughly 25 and 60% for protons and an- 
tiprotons respectively at pT — 0.5 GeV/c. 

Background-subtracted yields are corrected for angu- 
lar resolution of the DC and for smearing of the recon- 
structed momenta resulting from multiple scattering of 
tracks, which depends on hadron mass. To account for 
the acceptance of the PHENIX detector system and the 
varying efficiency for different hadrons, single-particle 
Monte-Carlo simulations are performed and verified by 
comparing the live detector area between data and Monte 
Carlo, and appropriate correction factors are determined 
separately for each hadron species. 

TABLE I: Fit-function parameters for the efficiency curves 
for different hadron species. See text for details. 



hadron ABC 

7r+ -0.08 -1.8 0.1453 

K+ -0.17 -0.97 0.1280 

p -0.21 -1.54 0.1428 

TT" -0.07 -1.7 0.1449 

R- -0.17 -1.01 0.1276 

p- -0.21 -1.56 0.1424 



Figure [2] shows the efficiencies for the three positive 
hadron species. The small efficiency for kaons is due to 
decays in flight. The large decrease in efficiency at low 
Pt is due to the fact that the fixed pseudorapidity accep- 



TABLE II: Fit-function parameters for relative fractions of 
different species in the hadron mix. See text for details. 



hadron 


A 


B 


C 


D 


7r+ 


1.02 


-2.39 


0.57 




K+ 


-0.53 


-2.39 


0.20 


0.009 


P 


-0.49 


-2.39 


0.23 


-0.009 




1.17 


-2.49 


0.61 


0.012 


!<- 


-0.61 


-2.49 


0.20 






-0.56 


-2.49 


0.18 


-0.012 



TABLE III: Systematic uncertainties of cross section mea- 
surements from various sources. 

Source Systematic Uncertainty 

Acceptance & efficiency correction 11-24% 

0"BBC 11% 

Trigger bias 2.5% 

Monte Carlo/data scale factor 2% 

PID fraction 1-5% 

Background fraction 1-5% 

Momentum smearing correction 1-2% 



tance of the detector corresponds to a narrow range in 
rapidity for smaller pT/m. The efficiencies are parame- 
terized (Ae^P^ -|-C) as a function of pt- TableUshows the 
fit-function parameters. The fit values are used for the 
calculation of cross sections. The species-dependent cor- 
rections are applied according to their production frac- 
tion (Fig. [21) in the hadron mixture. The production 
fractions were determined from identified hadron spec- 
tra from PHENIX [l^l, as well as earlier data from the 
ISR Dotted and dashed lines in Fig. [5] represent 

fits with individual data sets. Parameterized fits of the 
form Ae^P^ + C (for positive pious and negative kaons) 
and Ae^P'^ + C + DpT (for all other species) are per- 
formed under the constraint that the sum of the relative 
fractions is 1. Table HIl gives the fit- function parameters. 
Relative fractions from fit values are used to apply the 
corrections. The corrected yields are scaled by the BBC 
trigger bias as described in Ref. [ioj . 

The integrated luminosity Cint = "^f^i required for 
normalization of the invariant cross section, is calculated 
using the count of BBC-triggered events and the BBC 
normalization parameter cbbc- The parameter ctbbc is 
the p + p cross section seen by the BBC and is measured 
by the Van der Meer (Vernier) scan technique 0] . The 
quantity cbbc = 13.7 ± 1.5 ^^^^ mb for the relevant data 
set has been measured for previous PHENIX results and 
was discussed in detail in Ref. (40j . 

Table IIIII shows the systematic uncertainties of cross 
section measurements from various sources. The largest 
contribution (11-24%) to the pT-dependent systematic 
uncertainty comes from the correction for the acceptance 
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FIG. 3: (color online) Relative fraction of each species for positive hadrons. The error bars on the data points are combined 
statistical and systematic uncertainties. The solid line represents a fit to both the PHENIX and ISR data. The dashed (dotted) 
line represents a fit to only the PHENIX (ISR) data. 



and detection efficiencies. The uncertainty on the cross 
section due to this correction factor is determined by 
varying the selection parameters in the MC simulations. 
The trigger bias introduces a 2.5% uncertainty in the 
overall normalization, in addition to the 11% uncertainty 
on (Tbbc ■ Determination of the background fraction and 
the production fraction of separate hadron species each 
introduces a 1-5 % pr-dependent systematic uncertainty. 
Uncertainties from other sources, for example the correc- 
tion for momentum resolution, the correction for the ac- 
tive area of the detector in experiment and Monte Carlo, 
are ~l-2%. 



B. Cross Section Results 

Figure |3] and Table IIVI show the inclusive charged- 
hadron cross sections from p + p collisions at a/s = 62.4 
GeV as a function of pT- A combined pT-independent 
normalization uncertainty of 11.2% (uncertainties in the 
measurements of ctbbc ^^'^ BBC trigger bias) is not 
shown. 

In the overlapping px range, the results were found 
to be consistent with the species-combined cross sec- 
tions from identified results at PHENIX [TJ] as well as 
ISR results at ^/s = 63 GeV [i^. On the upper pan- 
els of both plots in Fig. |4] the cross sections are com- 
pared to NLO and NLL calculations at a factorization, 
renormalization and fragmentation scale of /i = pt [49| . 
The calculations were performed using Martin-Roberts- 
Stirling-Thorne (MRST2002) PDFs 50] and deFlorian- 
Sassot-Stratmann (DSS) fragmentation functions [sij . 
The NLO predictions have been shown to describe midra- 



pidity cross section results for neutral pious [52l, 'Hsj and 
charged hadrons [Hi at = 200 GeV within - 20% 
for a scale choice of = pt- For the present results 
at y/s = 62.4 GeV, however, the NLO calculations un- 
derpredict the data by as much as ^ 80% in the case 
of positive hadrons and ^ 60% in the case of negative 
hadrons for the same scale choice. The lower two panels 
in the Fig. 0] plots show the dependence of the theoretical 
calculations on the choice of factorization, renormaliza- 
tion and fragmentation scale (/i) for three different values 
{pT, Pt/2, and2pT)- The inclusion of higher-order terms 
in the NLL calculations leads to a considerably smaller 
scale dependence. 

These new data at an energy intermediate to typical 
fixed-target and collider energies are timely, as the details 
of how to work with resummation techniques in different 
kinematic regimes are currently being explored by the 
theoretical community (see for example [I,[l3,[i3|). Com- 
parison of the present results to the calculations at NLO 
with and without NLL terms included indicates that in 
the measured kinematic range, NLL terms make relevant 
contributions to the cross section. However, the tendency 
of the NLL calculations to overpredict the data, by as 
much as ~ 40% in the case of positive hadrons and ~ 50% 
in the case of negative hadrons for a scale choice of pr , 
may indicate that there are terms in the full next-to-next- 
to-leading-order (NNLO) expansion that are of compa- 
rable magnitude and opposite sign to those in the NLL 
calculation. These measurements corroborate similar in- 
dications from neutral pion cross section results [40| and 
identified hadron cross section results [l^ at PHENIX. 
The present measurements can also be useful in a future 
determination of inclusive charged-hadron fragmentation 
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FIG. 4: (color online) Cross section of inclusive-charged-hadron production at midrapidity in p + p at 62.4 GeV for (a) positive 
and (b) negative hadrons. NLO and NLL theoretical predictions [4^ at midrapidity, using MRST2002 parton distribution 
functions [50| and DSS fragmentation functions fsi!], at factorization, renormalization, and fragmentation scale jj. = px are 
shown as curves. The lower panels show the scale dependence of the NLO and NLL results. 



functions, as progress in pQCD has allowed inclusion of 
p+p cross section measurements and semi-inclusive deep- 
inelastic lepton-nucleon scattering data in FF parameter- 
izations alon g w ith the traditionally used e'^e" data since 

2007 [sil,!!!^. 



IV. DOUBLE-HELICITY ASYMMETRY 

We measured the double-heHcity asymmetries, All, of 
inclusive positive- and negative-hadron production in the 
transverse-momentum range of 0.5 < pt 4.5 GeV /c at 
midrapidity from longitudinally polarized p + p collisions 
at ^/i = 62.4 GeV. 



A. All Measurement 

The double-helicity asymmetry of charged hadrons 
is defined as the relative difference between hadron- 
production cross sections from collisions of the same- 
and opposite-helicity state protons. Experimentally, the 
asymmetry is measured as 



All = 



I N++ -R-N+. 

Pb ■ Py N++ +R-N+. 



(3) 



where Pb, Py are polarizations of the two colliding 
beams in RHIC (termed 'Blue' and 'Yellow'), 7V+_ 
are the midrapidity hadron yields from collisions of the 
same- and opposite-helicity protons, and relative lumi- 
nosity R — is the ratio of luminosity of the same- 
helicity collisions to that of opposite-helicity collisions. 

For the 2006 p + p data set at = 62.4 GeV, the 
luminosity-weighted average beam polarizations for both 
beams are measured to be (P) = 0.48, and the aver- 
age magnitude of the product of polarization of the two 
beams is {Pb ■ Py) — 0-23 with a relative uncertainty 
of 13.9%. The colliding proton bunches at RHIC are as- 
signed preset spin patterns repeated every four crossings. 
For consecutive fills with 120 bunches in the RHIC ring 
four such different spin patterns are alternated in order 
to reduce false asymmetries and systematic effects of pos- 
sible correlations between the detector response and the 
RHIC bunch structures. 

Hadron counts are obtained under similar criteria 
as described for the cross section measurements from 
the reconstructed tracks in Sect. IIII Al Approximately 
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TABLE IV: Cross section of midrapidity charged- hadron production from p + p collisions at ^/s — 62.4 GeV as a function of 
Pt- The errors represent the statistical (first) and systematic uncertainties. The data are corrected for the contribution of 
feed-down protons and antiprotons. A normalization uncertainty of 11.2% is not included. 
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1.63 X 10® BBC-triggered events with longitudinal beam 
polarization, corresponding to an integrated luminosity 
of 11.9 nb~^, were analyzed for the asymmetry measure- 
ments. Luminosities for the same- and opposite-helicity 
events were obtained from crossing-by-crossing informa- 
tion of BBC trigger counts. The systematic uncertainty 
of the relative luminosity i?, determined by comparing 
BBC-triggered events with events triggered by the ZDCs, 
was found to be 1.4 x 10~^. 



The asymmetry is determined on a fill-by-fill basis. 
The asymmetry of backgrounds from decays in flight, se- 
lected from the tail ends of the distribution of the match- 
ing variables (in beam direction z and azimuthal angle 0), 
is measured. The background asymmetry A^j^ and back- 
ground fraction are used to calculate the signal asymme- 
try and its statistical uncertainty. The feed-down from 
decays of A's and heavier hyperons, however, cannot be 
separated from hadron yields. The backgrounds from 
feed-down protons and antiprotons are a small contri- 
bution (1-2.5%) to the total hadron yields and are not 
corrected for since the asymmetries of these backgrounds 
are unknown. 



B. All Results 

Figure [S] and Table |V] show the p^-dependence of 
the measured double-helicity asymmetries for inclusive- 
charged-hadron production at midrapidity in polarized 
p+p collisions at ^/s = 62.4 GcV. The asymmetries 
are compared to NLO pQCD predictions [56] based on 
several different parameterizations of polarized PDFs at 
scale ^ = Pt- 

The curves in Fig. [S] labeled "DSSV NLO" and "DSSV 
NLL" refer to deFlorian-Sassot-Stratmann- Vogelsang 
(DSSV) parameterizations of the helicity PDFs [35| and 
curves labeled "BB NLO" refer to Bliimlein-Bottcher 
(BB) parameterizations of the helicity PDFs [sl]. Both 
DSSV and BB calculations use DSS fragmentation 
functions [H. DSSV calculations use MRST2002 
unpolarized PDFs 50]; BB calculations use those 
from the coordinated-theoretical-experimental-project- 
on-QCD-6 [13. Polarized PDFs use fits to the pDIS data 
to extract parameters of the functional forms of PDFs. 
DSSV [35] parameterizations use RHIC data along with 
the available pDIS data to constrain the polarized PDFs. 
The asymmetries are also compared to the NLL estima- 
tions with the DSSV PDFs [49]. 

For the purpose of comparison with the experimental 
results, pQCD calculations were obtained for separate 
hadron species (pions, kaons and (anti)protons) and were 
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FIG. 5: (color online) Double-helicity asymmetry (All) of (a) positive and (b) negative charged-hadron production from 
polarized p+p collisions at -^5 = 62.4 GeV. The results are compared to NLO and NLL pQCD predictions using several 
parameterizations of the helicity PDFs (see text for details). 



TABLE V: The double-helicity asymmetries and the statistical uncertainties as a function of pr for positive and negative 
inclusive charged hadrons from p + p collisions at -^s = 62.4 GeV. The fractional contribution to the yields from weak-decay 
feed-down to protons and antiprotons is shown; no correction to the asymmetries has been made for these contributions. 







estimated 




estimated 


PT (GeV/c) 




feed-down fraction(/i^) 


All + Mli,(/i") 


feed-down fraction(/i~) 


0.68 


0.0023 ± 0.0022 


0.022 


-0.0065 + 0.0024 


0.021 


1.26 


-0.01096 ± 0.0048 


0.016 


-0.0096 ± 0.0052 


0.025 


2.29 


0.0132 ± 0.0204 


0.012 


0.0208 + 0.0236 


0.021 


3.40 


-0.0299 + 0.0517 


0.011 


-0.0761 + 0.0578 


0.018 



combined using the particle fraction in the hadron mix- 
ture (Fig. [3]) and corresponding efficiency factor (Fig. [2]) . 
The measured asymmetries are small and consistent with 
zero. The results are also consistent with the predictions 
from the recent parameterizations within statistical limi- 
tations. The comparisons corroborate previous PHENIX 
measurements [4, 40] that disfavor very large gluon polar- 
ization. The presented asymmetry measurements probe 
a range of approximately 0.05 < Xgiuon < 0.2 [T^ of the 
interacting gluons. 



V. SUMMARY AND CONCLUSIONS 

Cross sections and double-helicity asymmetries for the 
midrapidity production of positive and negative inclusive 
charged hadrons at ^/s — 62.4 GeV are measured as a 
function of transverse momentum. The comparison with 
pQCD calculations shows that the NLO estimations are 
consistent with cross section results within a large scale 
uncertainty. NLL calculations with their reduced scale 
dependence are also consistent with the data, indicat- 
ing that the threshold resummation of logarithmic terms 
is relevant in the kinematic region measured; however, 



the overprediction of the data by up to ~ 50% if the 
NLL terms are included suggests that contributions from 
NNLO terms may also be important. This corroborates 
other recent results from PHENIX J^, i4Q,] with similar 
indications. The asymmetry results are the first measure- 
ments for charged hadron production in polarized p + p 
collisions at ^/s — 62.4 GeV and are consistent with the 
asymmetries found using several other probes at different 
collision energies at RHIC [1, d, i, [Ifl 0, HI ■ Experi- 
mental measurements of a variety of processes covering a 
broad kinematic range are essential to advancing our un- 
derstanding of QCD in hadronic interactions and nucleon 
structure, and the present measurements contribute to- 
wards that end. 
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